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INTRODUCTION 



Disc galaxies exhibit a variety of non-axisymmetric struc- 
ture (bars, spiral structure, lopsided structure, etc.). See 
Baldwin, Lynden-Bell & Sancisi (1980) and Richter & San- 
cisi (1994) for an HI survey of asymmetries in the discs of 
galaxies, and Rix & Zaritsky (1995) for a study of the phe- 
nomenon in stellar light. These asymmetries suggest the 
following general problem: what are the possible station- 
ary configurations of a two-dimensional self-gravitating fluid 
other than an axisymmetric razor-thin disc? One example 
is the family of Riemann discs (Weinberg & Tremaine 1983, 
Weinberg 1983) but these are uniformly rotating and hence 
only relevant to the central parts of disc galaxies. 

In this paper we address a modest component of this 
problem: we seek non-axisymmetric razor-thin discs of two- 
dimensional barotropic fluid that are stationary in an iner- 
tial frame. We distinguish between 'razor-thin' and 'two- 
dimensional,' applying the latter term to the equation of 
state (see Section 2). Furthermore we assume that our 
systems are scale-free, which reduces the partial differen- 
tial equations describing the system to ordinary differential 
equations. We also allow for the presence of an axisymmet- 
ric background potential. This simple and highly idealized 
problem already exhibits a rich variety of solutions. 

In Section 2 we derive the equations of motion. In Sec- 
tion 3 we solve these equations analytically for small de- 
partures from axisymmetry, and in Section 4 we construct 
highly non-axisymmetric discs numerically. 



ABSTRACT 

We develop equations to describe the equilibrium state of razor-thin scale-free 
barotropic fluid discs, with rotation curve v oc R~°, (3 £ (— \,\)- The discs may 
be embedded in a scale-free axisymmetric background potential. Nearly axisymmetric 
discs are constructed analytically and discs with azimuthal density variations as large 
as 10:1 are constructed numerically. For small departures from axisymmetry, we find 
that: (i) Stationary self-consistent cold or cool discs with m > 2 do not exist; (ii) on a 
plane whose axes are the Mach number and the strength of the background potential, 
there are generally two one-parameter sequences of non-axisymmetric discs (one with 
aligned isodensity contours and one with spiral contours); the spiral sequence exists 
only for m = 1; (iii) isolated discs with a flat rotation curve Q3 = 0) support non- 
axisymmetric equilibrium states at all Mach numbers, and discs with (3 = \ support 
a two-parameter family of self-similar spiral patterns. 

Key words: hydrodynamics - galaxies: kinematics and dynamics 



2 THE EQUATIONS OF MOTION 

We use 'razor-thin' to refer to a disc that has a density 

P (t) = Y.{R,0)S(z) (1) 

where (R, 6, z) are cylindrical polar co-ordinates, and S(z) 
is the Dirac 5-function. 

We use 'two-dimensional' to refer to a fluid in which 
the pressure II acts only in a plane. We use a barotropic 
equation of state of the form 



n = fcs 1 



(2) 



where k > 0, F > 0. 

The gravitational potential in the plane of the disc is 
given by the Poisson integral 



$(R,0) = -G 



E(r, <f>) rdr d</> 



[r 2 + R? - 2Rrcos(9 - <f>)] 



1/2' 



(3) 



If the disc is stationary, the velocity field v = vr&u + ve&e 
satisfies the continuity equation 



V ■ (Sv) = 

and the Euler equation 

(v V)v = - V($+¥ + #) 



(4) 



(5) 



where "3> represents an imposed axisymmetric background 
potential and H is the enthalpy. In a two-dimensional 
barotropic fluid the enthalpy is a function of the surface 
density, H = H(£). 
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Figure 1. (a) The locus of self-consistent solutions for (3 = 0.2, m = 1 (two leftmost solid curves) and m = 2 (right solid curve); 
the lower of the two m = 1 curves is the family of aligned discs and the upper is the family of spiral discs, which terminates at the 
bifurcation point (w, /) = (0.29432, 0.07910). The horizontal coordinate w is the inverse square Mach number and the vertical coordinate 
/ is the ratio of the fixed background potential to the unperturbed potential of the disc. For / < the background potential is repulsive 
and hence unphysical; the dashed line / = corresponds to an isolated disc with no background potential. In the shaded region, the 
axisymmctric disc is subject to an m = instability, (b) Same as (a), but for (3 = —0.1; the bifurcation point at which the m = 1 spiral 
family begins is (w,f) = (6.38318,0.04949). 
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Figure 2. The locus of non-axisymmetric equilibria for isolated 
discs, / = (the solid curve and the solid horizontal line at 
/3 = 0). In the shaded region the axisymmetric disc is subject to 
an m = instability. 

We now impose scale-invariance and reduce the two- 
dimensional partial differential equations to ordinary differ- 
ential equations. We write 

2 ^M, 



$ = -R 
where 

<p(R,6) = 8 + fi\ogR 



$ = -P~ 2/3 P, 



(6) 



(7) 



and (3 and \x are real constants. The disc has a spiral pattern 
if fj, 7^ 0; without loss of generality we can assume fi > 0. 
Equations (3), (4) and (5) can only be satisfied when 



v e = R- p b{<p), 



S = R~ 

v R = R^ f3 a(^), 
and 

H = R- 2f> Q{<p). 

Since we impose scale invariance, all physical quantities are 
determined by their values on the unit circle R = 1. Thus, 
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Figure 3. The location of the bifurcation point between the 
aligned and spiral families of m = 1 equilibrium solutions, 
as a function of j3. From left to right, the diamonds mark 
/3 = 0.4,0.3,0.2,0.1,0,-0.05,-0.1,-0.15. The curve is singu- 
lar (w — > oo) at f3 = —0.2071, where / = 0.4381; in a second 
branch with (3 — > — i the curve approaches w — > 0, / — > oo but 
lies beyond the boundary of this figure. 



for instance, we shall refer to E and S as 'surface density' 
interchangeably. 

Mass distributions with j3 > ^ are unphysical because 
they contain infinite point masses (J°° T,(R,9)RdRd8 — » oo 
as r — > 0); we summarise other constraints on /3 at the end 
of this section. 

The case /3 — is special and is most easily handled 
by taking limits as (3 -> 0, with P(ip) = P /{2/3) + Pi (</?), 
P = P /(2/3), Q(<p) ee Q /(2/3) + Qi(<p). The result is 



$ = 
S = 



Po log R — Pi (ip) + const, 
S(<P) 



<I> = Po log R + const, 



(8) 



R ' 

vr = a(ip), 
ve = b(ip), 

H = -Q logP + Qi(<p) + const. 
In practice we shall not need these equations, since all the 
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expressions for physical quantities below are well behaved in 
the limit (3 -> 0. 

The sound speed v s is given by 

2 dn 



^=rfE= fcrE " 
The enthalpy is given by 

To preserve scale-invariance we require that 
r= 1+4 ^. 



(9) 



(10) 



(11) 



1 + 2/3' 

the constraint T > implies (3 > — \. Equation (10) thus 
implies that 

Q(<p) - j^^- 1 = ^*sM 2/3/(2/3+1) , (12) 

Equations (3), (4) and (5) can be written 



P(<P) 



= G 

Jo [1 + 



y 2 - 2ycos((p - \ - P logy)] 1/2 ' 



3(3Sa - n{Sa)' - (Sb)' = 0, 



(13) 



(14) 



2/3(P + P-Q)-n(P- Q) - (3a -b z + fiaa + ab = 0,(15) 
and 

(P - Q)' + (/3 - l)ab - bb' - /j,ab' = 0, (16) 

where a dash denotes a derivative with respect to tp. Using 
the equation of state (12) to eliminate Q, these four equa- 
tions are to be solved for the four functions P, S, a and 
b. 

The allowable range of (3 is constrained by the con- 
vergence of Poisson's equation (3). Let us assume that 
S(tp) = s m ex.p(irmp) = s m R lm>1 exp(im9), m > 0. Then 
we show in the Appendix that 

Pbp) = Gs m Y m (0)e imv , -\m < (3 < ±(1 + m). (17) 
where 

= [3 - \imii, (18) 
and (Kalnajs 1971) 



Y m {0) 



r(f-0 + ±)r(f+0) 
r(f-0 + i)r(f+0 + i) 



(19) 



To determine the allowed range of (3 we first consider m/0. 
For m = 1 the allowed range of (3 is ( — |, 1) (outside this 
range the force from material at r — > oo or r — > diverges); 
for higher m the allowed range is broader. 

The case m = is special. The potential diverges and 
equation (19) does not apply unless (3 £ (0, |). However, it 
is easy to show that the total force arising from the surface 
density so is finite throughout the larger range (3 € (— \ ,h) 



and is given by 

F = -V$ = -20 R- 2 ?' Gs o Y a (0)e R , 



(20) 



the same as predicted by equation (17). In other words 
the range of validity of equation (17) can be extended to 
(3 € (—5, 5) for m = if it is only used to compute forces, 



as in equations (15) and (16). 

An additional constraint for warm discs (k > 0) is (3 > 
-j so that pressure increases with density (equations 2and 
11), which ensures that sound waves are stable. 

Henceforth we shall assume that (3 is restricted within 
the limits (—5, 5) for cold discs and (— |, 5) for warm discs. 



3 LINEAR THEORY 

We first examine the behaviour of scale-free discs that are 
slightly perturbed from axisymmetry. 



3.1 Axisymmetric state 

In an axisymmetric disc the variables are independent of the 
angle ip: 



P ^ = T0 P= T0 QM = W 

a(ip) = 0, b((p) = 6 , S(ip) = c. 



(21) 



The notation is chosen so that the radial force from the 
axisymmetric disc is —1 at R = 1; the parameter / > is 
the ratio of the force from the background potential to the 
force from the self-gravity of the disc. The surface density c 
is given by Poisson's equation (17): 

2f3GcY {/3) = 1. (22) 
The rotation speed bo is determined by equation (15): 



60 = 1 + / - 9- 



(23) 



The enthalpy is related to the surface density by the equa- 
tion of state (12): 



5 = (l + 4/3)fe C 2 ' 3/( ^ +1) . 



(24) 



We parametrise the temperature of the disc by the in- 
verse square of the Mach number (cf. equation 9), 



k 1 + 4/3 2/3/(2/3+1) 



b\R-^ b%\ + 2(3 



6§(l + 2/3)- 



(25) 



We shall characterise axisymmetric discs with a given 
rotation parameter (3 by the parameter pair (w,f), i.e. a 
measure of the temperature and a measure of the importance 
of the background potential. In terms of these parameters 
the rotation speed and enthalpy are given by 



1+/ 



(l + 2/3)(l + /> 



(26) 



l + (l + 2/3)w' J 1 + (l + 2/?)io 
where bg is given in terms of / and w by equation (26). 

3.2 Response to potential perturbation 

We impose a potential that is slightly perturbed from ax- 
isymmetry: 



(27) 



with p< 1, and m > an integer. We write the enthalpy 

as 



(28) 
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and the velocity as 

a(<p) = aie im ", b(i P ) = bo + b 1 e imlp 1 (29) 

from which equations (15) and (16) yield to first order 

20' (p — q) — 2&o&i + imboai = 0, 
im(p — q) + bo [(0 — l)oi — imbi] = 0, 



(30) 



where 0' is defined in equation (18). Equations (30) can be 
solved to find a\ and b\ in terms of p and q: 



ai 



2im(p - q)(0' - I) 
b [m 2 + 2(0 - 1)] ' 

(p-q) [m 2 + 2(3' '09-1)] 



(31) 



o [m 2 + 209 - 1)] 
The perturbed surface density may be written 

5 = c (l + o-e lmv ) (32) 

and to first order equation (14) yields 

(2/3' + /3)oi - im(&i + b <r) = 0. (33) 

From the equation of state (12) the perturbed enthalpy 
and surface density are related by 

1 +4/3 fcc 2/3/(2/3+l) 



? = 



1 + 2(3 

using equations (24) and (26) we find the simpler form 
q = b 2 wa. 



(34) 



(35) 

For a given unperturbed disc, specified by (3, w, and /, and 
a given azimuthal wavenumber m, equations (31), (33) and 
(35) determine the surface density and velocity perturba- 
tions induced by the imposed potential p. 

We now find the potential response to the perturbation, 
by applying equations (17) and (19) to equation (32): 



P r = Gc [Y o (0) + Y m (0')oe lmv ] 



(36) 



P r is the response potential to the imposed potential per- 
turbation Pi. Using equation (22) we may write equation 
(36) as 



Pr = ^+Ape^, 



(37) 



where we evaluate the Love number A using equations (22), 
(31), (33) and (35): 



Y m (0') 
2/% (/3) 



Abt 



m 2 + 2((3 - 1) 
m 2 + 2(3 + 2(3' - 4/3' 2 



(38) 



If A is positive then the response potential supports the per- 
turbing potential; if A is negative then the response poten- 
tial is anti-aligned with the perturbing potential, and cannot 
support it. The special case (3 = can be handled by taking 
the limit (3 — ► in equation (38) . 

3.3 Self-consistency 

If the perturbation supports the potential, then it may be 
possible to construct a self-consistent disc, in which P r = Pi. 
To do so we set A = 1, which requires 



Y m ((3') a 

— 0() 



2(3Y (/J) 



m 2 + 2((3 - 1) 
m 2 + 20 + 2(3' - 4(3' 2 



(39) 



To analyse the solutions of (39) we write the equation as 



B = bl(w-C) 1 



(40) 



where o 2 , is real but B and C may be complex. Thus we 
have 



%t(B) = bo[w - 5ft(C)], 3(B) = -b 2 Q(C). 



(41) 



The quantity Ss(C) vanishes only if p, = 0, or (3 = \ [or if 
m 2 — 2(1-/3), but this cannot occur for (3 6 ( — 5, 5)]. Thus 
we have three families of solutions:, (1) £s(C) = 0, p = 0; (2) 
9(C) / 0; (3) 9(C) = 0, (3 = \. We examine each family 
in turn. 



3.3.1 Aligned discs (p, = G). 

In this family the isodensity contours are aligned. The con- 
dition for self-consistency is derived by setting (3 = (3' in 
(39): 



Y m (J3) 
2(3Y ((3) 



= bo 



m 2 + 2((3-l) 
m 2 - 4(3((3 - 1) 



(42) 



For a given rotation parameter (3 and azimuthal wavenum- 
ber m, equation (42) provides a unique relation between 
the Mach number (parametrised by to) and the back- 
ground potential (parametrised by /) that support a non- 
axisymmetric stationary state; in other words the non- 
axisymmetric aligned discs are located on a curve in the 
(w, f) plane. Figures 1 and 2 plot solutions of equation (39) 
for various values of (3 and m. 

The qualitative characteristics of these solutions are il- 
luminated by examining special cases: 

(i) Cold discs (w = 0). The left side of (42) is positive 
for all (3, so self-consistency requires 



m 2 + 2{(3 - 1) 
m 2 - 4(3(0 - 1) 



<0, 



(43) 



which is satisfied for m > and £ (— i, i) if and only if 
m = 1 and > -±(2 1/2 - 1) = -0.2071. A more stringent 
condition is that / > since a repulsive background is un- 
physical; for m = 1 and w = 0, equation (42) yields / > 
only when > 0. Thus cold, non-axisymmetric stationary 
states can only exist if m = 1 and > 0. 

(ii) Low-mass discs (/ — > 00). In this case the quan- 
tity in square brackets in equation (42) must be zero, which 
requires 



m 2 + 2Q3 - 1) _ 
rn 2 - 40(0 - 1) ' 



(44) 



since w > 0, stationary non-axisymmetric discs with neg- 
ligible mass exist if and only if the inequality (43) is not 
satisfied; in other words if and only if m > 2 or < 
-i(2 1/2 - 1) = -0.2071. A corollary is that for m = 1 
and > —0.2071, / is always finite; i.e. the disc always has 
significant self-gravity. 

(iii) Non-rotating discs (w — ► 00). In this case 



(l + 20)Y m (0) 
1 20Yo(0) 



(45) 



For m = 1 the condition / > requires < 0. For m > 2 
equation (45) yields unphysical solutions (/ < 0) for all 0. 
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(iv) Flat rotation curve (/3 = 0) . We may evaluate (42) 
using the relation T(x) ~ x^ 1 as x — > 0; in particular for 
m = 1 we find the simple result that 



f — for /3 = 0, and all w; 



(46) 



in other words isolated discs with a flat rotation curve sup- 
port a neutral m = 1 mode for any Mach number. 

Comparison of (iii) with (i) shows that the in — 1 so- 
lution curve in the (w, f) plane must cross / = for all 
rotation parameters /3; the location of this crossing point 
(i.e. the value of w that supports isolated non-axisymmetric 
discs) is shown in Figure 2. 



3.3.2 Spiral discs (fj, > 0) 

These are discs with self-similar spiral patterns. In equation 
(41), 3(C) 7^ 0. Equation (39) has real and imaginary parts. 
Thus, once (3 and m are specified, there are two equations 
for three unknowns (w, f,fi). Therefore the spiral discs are 
located on a curve in the (w, f) plane parametrised by the 
value of fi. The aligned and the spiral sequences are plotted 
in Figure 1. 

The spiral sequence bifurcates from the sequence of 
aligned discs at the point where fi — > 0. The location of 
the bifurcation point can be determined analytically. For 
small /i we can replace Q(B) by — ^irnii(dB/df3') where the 
derivative is evaluated at (i = 0, with a similar expression 
for 9(C). Thus at the bifurcation point, the equilibrium 
equations (41) simplify to 



r .2, ~ dB 2 8C 
B = b (w-C), W = ~bo W , 



(47) 



where all quantities are evaluated at fi = (/3 = (3'). Elimi- 
nating bo, we find 



HogB _ <91og(™ - C) 



d(3' 



d(3' 



(48) 



the derivatives are taken with (3, f and w fixed. Evaluating 
this expression yields 



2(l-4/3)[m 2 + 2(/3-l)] m 2 + 2(/3 - 1) 



[m 2 + 4/3(1 - /3)] 2 tf m (/3) m2+ 4/3(1 - /3) : 
where 

* m (z) = ij){\m + z) + 4>(\rn - z + 1) 

— y>(§m + ;z+i)- V"(| m — z + 



(49) 



(50) 



and ip(z) is the digamma function. This expression defines w 
at the bifurcation point; the other coordinate / then follows 
from (42). The expression (49) is singular for m = 1 at 
(3 = -±(2 1/2 - 1) = -0.2071, / = 0.4381. 

The location of the bifurcation point is plotted in Figure 
3 for m = 1. For m > 2 the bifurcation point is unphysical 
(/ < 0) for all (3. 

The equation governing self-similar spiral patterns sim- 
plifies in the WKB limit (fi 3> 1). Using the asymptotic 
properties of the Gamma function we may show that in this 
limit 



Ym(P') = F m (/3+Im M ) = — f 1 + ^-^ 
1 m/j, \ imfi 



+0(/x- 3 );(51) 



thus the dominant real and imaginary parts of equation (39) 
become respectively 

m 2 + 2(13-1)' 



TV 
7T 



2bl 



V + 2(/3 - 1)]. 



(3Y (f3) mil- 
Taking the ratio of these two equations yields 
2(1 - f3) -m 2 



as [i 



Similarly, 
/ 



2(3Y (P) 



m/j, 



2(1-/3) 



as /i 



(52) 



(53) 



(54) 



Thus w — > and / — > oo as fi — > oo; in other words the 
WKB limit applies only if discs are cold and the background 
dominates the unperturbed potential. Moreover the WKB 
solution is unphysical for m > 2 since these equations imply 
that w and / are negative. Equations (53) and (54) can also 
be combined to eliminate fi: 



Wf = 



2/3Fo(/3) 



2(1-/3) -m 2 ' 



/ > 1.(55) 



It is informative to relate the properties of these discs 
to the usual WKB spiral density waves in a fluid disc (e.g. 
Goldreich & Tremaine 1979). The dispersion relation for 
these waves has two branches ("long" and "short" waves) 
so a typical disc can support zero, one, or two wavetrains 
at any point. Why then do self-similar spiral patterns only 
exist in discs satisfying the constraint (55)? Consider a self- 
similar disc with unperturbed surface density oc i?~ 2,3 ~ 1 . In 
this disc, linear WKB density waves with zero pattern speed 
are also self-similar, but the surface density perturbation as- 
sociated with the waves is generally oc B7 Z ^ , a result which 
follows from conservation of angular momentum flux (e.g. 
Goldreich & Tremaine 1979). Thus the fractional density 
perturbation is oc R~ 2l3+1 / 2 , which is not independent of ra- 
dius, so the combination of axisymmetric disc and density 
wave is not self-similar. The only exceptions are (i) (3 = j, 
a case which we examine below; (ii) waves with zero angular 
momentum flux, the condition for the existence of which can 
be shown to be the same as (55). 

We have seen that both the start of the spiral sequence 
(/i = + ) and the asymptotic limit of the sequence (n — ► oo) 
lie in unphysical regions of the (w, f) plane (w < or / < 0) 
when m > 2. These results suggest that the entire spiral 
sequence is unphysical when m > 2. We have not been 
able to prove this conjecture, but a numerical search has 
not revealed any counter-examples. 

3.3.3 S3=\ discs 

The condition for self-consistency can be written 

r(i m + i + i; m/i )r(f) 2 



r(im + | + iim M )r(i) 

1 + / 



2 + 3w 



w + ■ 



f-m 2 



(56) 



For this particular value of (3, equation (39) is real for all /it. 
Once m is specified, there is one equation for three unknowns 
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Figure 4. As in Figure 1, for a disc with (5 = -g. The solid curves 
show the loci of solutions with m = 1; from top to bottom at large 
w: n = 0, 1, 2, 4, 8, 16. The dot-dashed curve shows m = 2 and 
/i = 0. For m = 2 and fi > spiral solutions are unphysical 
because ui < 0; the region below the dashed line at / = is also 
unphysical. In the shaded region the axisymmetric disc is subject 
to an m = instability. 



(to, /, n). Thus the spiral discs are located on a set of curves 
in the (to, /) plane labelled by the value of /i. Self-similar 
spiral discs with m = 1 are not restricted to a single curve 
in the (to, f) plane: instead they occupy a two-dimensional 
area, as shown in Figure 4. 



3.4 Stability 

The linear stability of axisymmetric scale-free discs to ax- 
isymmetric disturbances has been investigated by Lemos 
etal. (1991). Their analysis is based on the observation that 
the squared eigenfrequency of axisymmetric modes is real, 
so that neutral modes mark the boundary between stable 
and unstable modes (Lynden-Bell & Ostriker 1967). Thus 
their calculations are quite similar to ours, except that they 
focus on the azimuthal wavenumber m = 0, whereas our 
analysis beyond equation (36) is only valid for m > 0. We 
may extend our analysis to the axisymmetric case by replac- 
ing exp(imtp) in equations (27), (28), (29), (32) and (36) by 
exp( if log R) and carrying out a slightly more careful deriva- 
tion using Lagrangian perturbations. In our notation, the 
condition for a neutral, self-similar, axisymmetric distur- 
bance is (cf. equation 39) 



W) 

2(3Y (/3) 



i 2 



/3' - 2/3' 2 



(57) 



where (3' = f3 — \iv. Since we are only investigating stabil- 
ity rather than constructing a self-similar non-axisymmetric 
disk, the fractional amplitude of the perturbation need not 
be independent of radius; thus f3' can be an arbitrary com- 
plex number with < K(/3') < ^ (equation 17). Following 
Lemos etal. (1991), we focus on 5R(/3') = |, for which both 
sides of equation (57) are real, although roots with other 
values of 5K(/3') may exist. 

The location of the stability boundaries in the (to, /) 
plane can be determined by the following argument. Let 
/3' = \ — \icx, where a is real, and rewrite (57) in the form 



where 
B{a) = 

and 

(7(a) = 



^|r(|-^a)| 2 

2/3Y (/3)|r(! + iza)| 2: 



2(/3 - 1) 



(59) 



(60) 



Since F(a) is even, it must have an even number of roots. 
The stability boundary is marked by the transition from zero 
roots to two roots, which occurs when 



F(a) = B(a) - bl[w - C{a)] = 0, 
F'(a) = B' (a) + b 2 C' (a) =0 

simultaneously for some a. Eliminating bo, we have 

dlogB(a) _ dlog[to - C(q)] 
da 

which yields 
i*o(| - \ia) = 



da 



8(1 - f3)a 



(61) 



(62) 



(63) 



(i + Q 2)[(i+a2) w + 2(1 -/?)]' 

where ^ m (z) is defined in equation (50) and both sides are 
real. To locate the stability boundaries we distinguish two 
cases: (i) if a = equation (63) is satisfied, and the corre- 
sponding stability boundary on the (w, /) plane is found by 
solving F(0) = for / as a function of w; (ii) if a =fc we 
can easily manipulate (63) to solve for w(a) and then use 
F(a) = to find / as a function of a and io(a). Figures 
1, 2, and 4 show the regions where m = instabilities are 
present in the (w,f) plane. 

The close relation between neutral modes and the on- 
set of axisymmetric instability does not extend to non- 
axisymmetric disturbances (e.g. Lynden-Bell & Ostriker 
1967). However, non-axisymmetric neutral modes would 
mark the stability boundary for systems composed of two 
equal, counter-rotating discs (we imagine that the discs 
interact only through their mutual gravity), since in this 
case the equilibrium is time-reversal invariant. Thus anal- 
yses similar to that of Lemos et al. could illuminate non- 
axisymmetric two-stream instabilities in scale-free counter- 
rotating discs (Araki 1987, Sellwood & Merritt 1994). 



4 NONLINEAR THEORY 

We may extend our analysis to nonlinear departures from 
axisymmetry by expanding all the disc properties in a 
Fourier series, truncating the series, and solving the re- 
sulting nonlinear equations for the Fourier coefficients by 
Newton's method. For a given rotation parameter /3 and 
background potential /, we expect to have a sequence of 
solutions with m-fold symmetry in which the Mach number 
parameter to is a function of the amplitude 



S(tt)/S(0) 



(64) 



where zero azimuth is chosen so that Z > 0. In the limit 
Z — > 0, w is determined by the linear theory of Section 3. 

We have restricted our numerical calculations to aligned 
discs. Thus we set 



F(a) = B(a) - b 2 [w - C(a)] = 



(58) 



S(0) = c 



1 + Sfe cos(fcm#) 



(65) 
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Figure 5. Cold non-axisymmctric discs with fi = 0.33. a) and b) show streamlines for density contrast Z = 2.333 and Z = 9 respectively, 
(c) and (d) show density contours for those discs. 



a(8) — ^""^ afc sin(km8), b(6) = &o + ^^ frfc cos(fcmfl)(66) 



fc=i 



fe=i 



and truncate all the series at order n. The potential of the 
disc is given by 



P{9) = Gc 



YoW) + Y km{/3)s k cos(km8) 



and the enthalpy is 

n 

Q(8) = ^+^2<lk cos(kme). 



(67) 



(68) 



For computational convenience we work with warm discs 
whose rotation parameters have the form /3= \{i — 
where £ is an integer (e.g. (3 = |, ^, |, . . .) since (12) then 
implies that S oc Q e , which is easier to evaluate if I is an 
integer. The solutions are parametrised by the value of Z. 



Figure 5 shows streamlines and density contours for two 
cold discs with j3 = 0.33, using m = 1, n = 5. The first 
has Z — 2.333 and / = 0.658. The second has Z = 9 and 
/ = 0.816. For comparison, linear theory predicts / = 0.604. 

Figure 6 shows streamlines and density contours for two 
isolated (/ = 0) warm discs with /3 = |, using m = 1 and 
n — 10. The first has Z = 2.333 and w = 0.355. The second 
has Z — 9 and w — 0.349 (linear theory predicts w = 0.353). 

Figure 7 shows streamlines and density contours for two 
isolated (/ = 0) cold discs with /3 — 0, using in — 1 and 
n = 10. The first has Z = 1; the second has Z = 2.333. 

5 DISCUSSION 

We have sought non-axisymmetric equilibrium states of 
scale-free razor-thin discs of a two-dimensional barotropic 
fluid embedded in an axisymmetric background potential. 
The rotation curves of the discs have the form v oc R^ 13 , 
G (—5,5) for cold discs and /3 £ (— 1,|) for warm 
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Figure 6. Warm isolated discs (/ = 0) with /3 = A. a) and b) sh 
density contours for those discs. 

barotropic discs. The discs are specified by the rotation 
parameter P, the inverse-square Mach number w > 0, and 
the ratio of the force from the background to the force from 
the unperturbed disc, / > 0. For small departures from 
axisymmetry, we have found that: 

(i) In cold or cool discs (w <C 1), stationary equilibrium 
states exist only for lopsided modes (m = 1) and for 
rotation parameters p > -§(2 1/2 - 1) = -0.2071. 

(ii) For most values of /3, stationary equilibrium states ex- 
ist only along two sequences in the (to, /) plane: the 
aligned discs and the spiral discs. The spiral-disc se- 
quence bifurcates from the aligned-disc sequence. We 
have argued that the spiral-disc sequence is only phys- 
ical if m = 1. 

(iii) Isolated (/ = 0) discs with flat rotation curves (J3 — 
0) support non-axisymmetric equilibrium states for all 
values of the Mach number. 

(iv) Discs with p — \ support m = 1 spiral equilibrium 
states not just along a single sequence but over a large 




-1 -0.5 0.5 1 







-0.5 




-1 -0.5 0.5 1 

streamlines for Z = 2.333 and Z = 9 respectively, (c) and (d) show 



part of the (w, f) plane. 
We have also constructed equilibria with large departures 
from axisymmetry (range of surface density at a given radius 
of a factor ten). 

The rotation curves of our models have the character- 
istic feature that they are asymmetric along the long axis. 
The rotation velocity is high where the surface density is low. 
Similar features are seen in the the rotation curves of Mag- 
ellanic irregulars (de Vaucouleurs & Freeman 1972). Other 
lop-sided galaxies, such as NGC 891, where the gas distri- 
bution extends much further on one side than the other, 
have symmetric rotation curves (Sancisi & Allen, 1979). In 
NGC 891 the stellar distribution is also very symmetric. Our 
models do not reproduce this behaviour. 

It would be interesting to construct stellar dynamic 
equivalents of our fluid models. There is no general argu- 
ment which guarantees that such models exist. However, 
the cold fluid discs we have constructed are also valid mod- 
els for cold stellar discs; also, the case P = 0, m = 2 has been 



Non-axisymmetric discs 9 





Figure 7. Cold isolated discs with fi = 0. a) and b) show streamlines for Z = 1 and Z = 2.333 respectively, (c) and (d) show density 
contours for those discs. 



studied by Kuijken (1993), and a number of hot models were 
constructed. 
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APPENDIX 

We wish to derive the potential arising from a surface density 
of the form 

E(i?,0) = R a e lm0 . (69) 
We begin with the potential-density pair (Qian 1992) 
(m+ i) (R\ m e im9 



ib 2 

G (R^ 



(1 + ^2/52)^+3/2 ' 



*.<*')=-!(!) 



(70) 



(1 + #762)m+l/2 • 
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We now form the scale-free density-potential pair 

b a+1 E b (R,9)db; / b a+1 <S> b (R,9)db; (71) 
Jo 



both integrals converge so long as — m — 2 < 5R(a) < m — 1, 
and yield 

1 „c i^r(im+ia + |)r(im-ia) 
—-H e 



2tt T(m+i) 

2 (72) 
_G +1 <m9 r(lm+|a + l)r(lm-la-l) 

2 K r(m+|) 

Thus the potential corresponding to the surface density (69) 
is (Kalnajs f97f) 

$(ii,6>) = -TvGR a+1 e tmt> 

r(l m+ l a + l)r(lm-l Q -l) (73) 
r(I m +I a +|)r(im-iQ) ' 

where — m — 2 < 5R(a) < m — 1. 
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